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ABSTRACT 

Type la supernovae have been proposed to be much better distance indicators at near- 
infrared compared to optical wavelengths - the effect of dust extinction is expected 
to be lower and it has been shown that SNe la behave more like ‘standard candles’ at 
NIR wavelengths. To better understand the physical processes behind this increased 
uniformity, we have studied the Y, J and H-Slter light curves of 91 SNe la from 
the literature. We show that the phases and luminosities of the first maximum in 
the NIR light curves are extremely uniform for our sample. The phase of the second 
maximum, the late-phase NIR luminosity and the optical light curve shape are found to 
be strongly correlated, in particular more luminous SNe la reach the second maximum 
in the NIR filters at a later phase compared to fainter objects. We also find a strong 
correlation between the phase of the second maximum and the epoch at which the SN 
enters the Lira law phase in its optical colour curve (epochs ~ 15 to 30 days after B 
band maximum). The decline rate after the second maximum is very uniform in all 
NIR filters. We suggest that these observational parameters are linked to the nickel and 
iron mass in the explosion, providing evidence that the amount of nickel synthesised 
in the explosion is the dominating factor shaping the optical and NIR appearance of 
SNe la. 
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1 INTRODUCTION 

The uniformity of Type la supernovae (SNe la) has led to 
their use, after calibration, as distance indicators (reviewed 
in: Goobar & Leibundgut 2011) and they provided the first 
evidence for the accelerated expansion of the Universe (Riess 
et al. 1998; Perlmutter et al. 1999). 

Observations of large SN la samples show that the peak 
luminosity in the optical is not uniform but can be nor¬ 
malised following a variety of calibration techniques, most 
notably the correlation between light curve shape and peak 
luminosity, and between light curve colour and peak lumi¬ 
nosity (e.g. Phillips 1993; Riess, Press, & Kirshner 1996; 
Guy et al. 2005, 2007, 2010; Jha, Riess, & Kirshner 2007). 
The variation in bolometric luminosity for the objects (Gon- 
tardo, Leibundgut & Vacca 2000) implies variations in the 
physical parameters of the explosion, in particular the syn¬ 
thesised Ni mass and the total ejected mass (Stritzinger et 
al. 2006; Scalzo et al. 2014). 
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At near-infrared (NIR) wavelengths (900nm > A > 
2000nm), SNe la have a very uniform brightness distri¬ 
bution without any prior normalisation (Elias et al. 1981; 
Meikle 2000; Krisciunas et al. 2004, 2007). The scatter in the 
peak luminosity in these studies is ~0.2 mag, which when 
combined with the lower sensitivity of the NIR to extinc¬ 
tion by dust, has sparked interest in the use of this wave¬ 
length region. Following large observational campaigns, sta¬ 
tistically significant samples of SN la light curves have been 
made public (Wood-Vasey et al. 2008; Gontreras et al. 2010; 
Stritzinger et al. 2011; Barone-Nugent et al. 2012), and have 
been used to construct the first rest-frame NIR Hubble dia¬ 
grams (Nobili et al. 2005; Freedman et al. 2009; Kattner et 
al. 2012; Weyant et al. 2014). 

Dust extinction in the NIR is significantly reduced com¬ 
pared to the optical leading to smaller corrections and un¬ 
certainties. In addition to the photometric calibration sys- 
tematics (see Conley et al. 2011) the dust extinction for SNe 
la is one of the major sources of systematics in SN la cosmol¬ 
ogy measurements (e.g. Peacock et al. 2006; Goobar & Lei¬ 
bundgut 2011). In particular, extinction law measurements 
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for galaxies remains uncertain (see discussions in Phillips et 
al. 2013; Scolnic et al. 2014). The recent SN 2014J is a case 
in point with derived dust properties very different from 
local Interstellar dust (Amanullah et al. 2014; Foley et al. 
2014). Strongly reddened SNe la may also exhibit variations 
in their light curve shapes (Leibundgut 1988; Amanullah et 
al. 2014). 

The light-curve morphology in the NIR is markedly dif¬ 
ferent from that in the optical with a pronounced second 
maximum in lYJHK filters for ‘normal’ SNe la (Elias et 
al. 1981, 1985; Leibundgut 1988, 2000; Meikle 2000; Wood- 
Vasey et al. 2008; Folatelli et al. 2010). The formation of the 
NIR spectrum in SNe la is highly sensitive to the opacity 
variations (Spyromilio, Pinto, & Eastman 1994; Wheeler et 
al. 1998) as the spectrum is dominated by line blanketing 
opacity making the evolution of the NIR light a sensitive 
probe of the structure of the ejecta. Kasen (2006) in a de¬ 
tailed study suggested that the second maximum is a result 
of a decrease in opacity due to the ionisation change of Ee 
group elements from doubly to singly ionised atoms, which 
preferentially radiate the energy at near-IR wavelengths. A 
direct prediction of Kasen (2006) is that a larger iron mass 
leads to a later NIR second maximum. 

Studies of the i band light curve find a relation between 
the phase of the second maximum and the optical light curve 
shape (e.g. Ami 5 (R), Folatelli et al. 2010; Hamuy et al. 
1996). The strength of the second maximum in i does not 
show such a correlation. 

In this paper, we investigate the properties of SN la NIR 
light curves (Y JH) and establish correlations with other ob¬ 
servational characteristics. Connections to possible physical 
properties in the explosions are explored. The structure of 
this paper is as follows: after a presentation of the input data 
in Section 2, we analyse the NIR light curve properties (Sec¬ 
tion 3) along with a description of NIR colours. Correlations 
with optical light curve parameters and their interpretations 
are given in §4 followed by a discussion in §5. The conclu¬ 
sions are presented in Section 6. 


2 DATA 

We investigate a large sample of nearby objects with well- 
sampled optical and NIR data (Table 1). The main data 
source of NIR SN la photometry is the Carnegie SN Project 
(CSP; Contreras et al. 2010; Burns et al. 2011; Stritzinger et 
al. 2011; Phillips 2012; Burns et al. 2014). The low-redshlft 
CSP provides a sample of SNe la with optical and NIR light 
curves in a homogeneous and well-defined photometric sys¬ 
tem (in Vega magnitude system) and thus forms an ideal 
basis for the evaluation of light curve properties. CSP re¬ 
lies primarily on SN discoveries from the Lick Observatory 
SN Search (LOSS; Leaman et al. 2011). The CSP has pub¬ 
lished light curves on a total of 82 SNe la of which 70 have 
photometry in YJHK bands. 

From the CSP NIR dataset, we removed spectroscopi¬ 
cally peculiar objects such as SN2006bt and SN2006ot. We 
also rejected SNe la with spectra similar to the peculiar 
SN 1991bg (Filippenko et al. 1992; Leibundgut et al. 1993; 
Mazzali et al. 1997) and objects that do not exhibit a second 
maximum (SNe 2005bl, 2005ke, 2005ku, 2006bd, 2006mr, 
2007N, 2007ax, 2007ba, 2009F). 


Table 1. SN sample 


SN name 

Phase range 
(days) 

Ni 

Ny 


Nfj 

"k 

Reference 

SN1980N 

5.7 

• • • 99.7 



11 

13 


MOO 

SN1981B 

2.7 

- • . 120.5 



17 

17 


MOO 

SN1986G 

-6.1 

• • • 101.0 



28 

29 


MOO 

SN1998bu 

-8.5 

- • • 31.5 



23 

23 


MOO 

SN1999ac 

-13.1 

• • • 63.8 



30 

30 


Ph06 

SN1999ee 

-7.5 

• • • 27.5 


17 

18 

20 


K04a 

SN1999ek 

-8.0 

- • • 23.0 



14 

15 


K04a 

SN2000E 

-7.7 

• • • 126.6 



18 

18 


V03 

SN2000bh 

-5.5 

- • • 39.0 


6 

21 

22 


K04a 

SN2001ba 

-6.0 

• • • 34.9 



14 

15 


K04a 

SN2001bt 

-1.6 

- • • 65.3 



21 

21 


K04a 

SN2001cn 

4.2 

• • • 59.1 



19 

19 


K04b 

SN2001CZ 

-2.4 

- • . 47.6 



12 

12 


K04b 

SN2001el 

-10.6 

64.3 



33 

32 


K03 

SN2002bo 

-11.0 

. . . 44,0 



17 

17 

17 

K04b, ESC 

SN2002dj 

-11.0 

- • • 67.0 



21 

21 


P08, ESC 

SN2002fk 

-12.2 

• • • 102.1 



24 

23 

23 

C13 

SN2003cg 

-6.4 

- • . 413.5 



13 

13 


ER06, ESC 

SN2003hv 

1.2 

• • • 62.0 


16 

16 

16 


L09, ESC 

SN2004ef 

-8.7 

- • . 65.2 

46 

4 

3 

4 

3 

CSP 

SN2004eo 

-12.0 

• • • 63.0 

39 

8 

9 

9 

8 

CSP, POT, ESC 

SN2004ey 

-8.9 

- • • 48.1 

32 

7 

9 

9 

8 

CSP 

SN2004gs 

-3.6 

• • • 99.1 

50 

12 

11 

10 


CSP 

SN2004gu 

-0.4 

- • . 48.6 

27 

8 

7 

7 


CSP 

SN2005A 

-3.5 

- • • 30.4 

35 

10 

10 

10 


CSP 

SN2005M 

-8.0 

- • • 74.7 

56 

17 

17 

14 

12 

CSP 

SN2005ag 

-1.7 

- • • 66.3 

44 

9 

9 

9 


CSP 

SN2005al 

-1.0 

• • • 82.0 

35 

7 

8 

8 

7 

CSP 

SN2005am 

-4.6 

- • • 75.3 

36 

6 

6 

6 

6 

CSP 

SN2005el 

-7.3 

• • • 83.6 

25 

21 

22 

15 

3 

CSP 

SN2005eq 

-3.6 

- • • 97.3 

27 

15 

15 

10 

1 

CSP 

SN2005hc 

-4.5 

• • • 84.4 

22 

13 

11 

9 


CSP 

SN2005hj 

-2.3 

- • • 87.6 

16 

11 

12 

10 

1 

CSP 

SN2005iq 

-5.3 

• • • 67.7 

19 

11 

11 

11 

1 

CSP 

SN2005kc 

-10.4 

- • • 30.6 

13 

9 

9 

8 


CSP 

SN2005ki 

-9.9 

• • • 155.8 

47 

12 

11 

10 


CSP 

SN2005na 

-1.8 

- • . 94.0 

27 

14 

11 

12 


CSP 

SN2006D 

-6.0 

- • • 112.8 

42 

17 

16 

16 

6 

CSP 

SN2006X 

-11.0 

• • • 119.8 

39 

32 

33 

32 

9 

CSP, ESC 

SN2006ax 

-11.7 

- • • 71.2 

26 

19 

18 

16 

3 

CSP 

SN2006bh 

-4.9 

• • • 60.0 

24 

12 

11 

10 


CSP 

SN2006br 

5.8 

- • • 37.7 

9 

5 

5 

5 


CSP 

SN2006ej 

4.4 

- • • 69.4 

13 

3 

3 

3 


CSP 

SN2006eq 

0.2 

- • . 43.0 

18 

10 

7 

8 


CSP 

SN2006et 

-7.2 

• • • 105.7 

23 

18 

12 

13 


CSP 

SN2006ev 

4.2 

- • • 48.1 

12 

10 

8 

8 


CSP 

SN2006gj 

-1.8 

• • • 96.0 

19 

13 

10 

4 


CSP 

SN2006gt 

-2.1 

• • • 62.8 

13 

10 

8 

6 


CSP 

SN2006hb 

7.2 

- • • 136.9 

25 

10 

10 

9 


CSP 

SN2006hx 

-9.3 

• • • 30.7 

8 

7 

6 

5 


CSP 

SN2006is 

6.4 

- • • 118.2 

24 

8 

8 

7 


CSP 

SN2006kf 

-4.8 

• • • 86.0 

20 

17 

14 

11 

3 

CSP 

SN20061U 

5.1 

- • . 90.1 

21 

6 

4 

3 


CSP 

SN2006ob 

-2.6 

• • • 58.3 

13 

12 

9 

5 


CSP 

SN2006OS 

2.4 

- • • 57.3 

14 

10 

6 

5 


CSP 

SN2007A 

-4.7 

• • • 18.3 

9 

9 

5 

3 


CSP 

SN2007S 

-10.5 

- • • 102.2 

19 

12 

17 

18 

7 

CSP 

SN2007af 

-8.8 

• • • 82.1 

28 

26 

25 

24 

5 

CSP 

SN2007ai 

-0.4 

- • • 86.4 

17 

7 

7 

6 

3 

CSP 

SN2007as 

-0.0 

- • • 77.7 

19 

11 

10 

10 


CSP 

SN2007bc 

-3.0 

• • • 58.8 

10 

11 

8 

6 


CSP 

SN2007bd 

-9.5 

- • . 215.8 

14 

12 

9 

7 


CSP 

SN2007bm 

-0.2 

• • • 34.7 

10 

10 

9 

7 


CSP 

SN2007ca 

-8.8 

- • • 32.1 

12 

10 

8 

7 

2 

CSP 

SN2007if 

14.4 

- • • 100.2 

16 

8 

7 

5 


CSP 

SN2007jg 

-3.1 

- • • 58.8 

18 

8 

6 

5 


CSP 

SN20071e 

-11.5 

71.2 

26 

17 

17 

16 


CSP 

SN2007nq 

-1.4 

- • • 77.5 

25 

19 

10 

4 


CSP 

SN2007on 

-8.6 

• • • 88.2 

38 

29 

28 

25 

7 

CSP 

SN2008C 

2.9 

• • • 85.6 

19 

15 

13 

18 

1 

CSP 

SN2008R 

-1.3 

- • • 32.6 

12 

8 

7 

5 


CSP 

SN2008bc 

-9.3 

• • • 96.6 

32 

7 

12 

11 


CSP 

SN2008bq 

-0.8 

- • . 43.1 

16 

4 

4 

4 


CSP 

SN2008fp 

-6.0 

• • • 89.9 

28 

22 

20 

20 

7 

CSP 

SN2008gp 

-10.5 

- • • 34.4 

19 

10 

11 

9 


CSP 

SN2008hv 

-10.5 

- • • 78.3 

25 

18 

16 

16 


CSP 

SN2008ia 

-3.5 

- • . 43.5 

15 

16 

15 

14 


CSP 

PTF09dlc 

-5.56 

• • • 19.5 



4 

4 


BN12 

PTFlOhdv 

-3.9 

- • • 10.3 



4 

4 


BN12 

PTFlOhmv 

-10.1 

7.9 




5 


BN12 

PTFlOmwb 

-12.3 

5.7 



5 

5 


BN12 

PTFlOndc 

-3.9 

5.1 



4 

4 


BN12 

PTFlOnlg 

-5.1 

3.9 



3 

5 


BN12 

PTFlOqyx 

-3.5 

7.5 



4 

4 


BN12 

PTFlOtce 

-6.6 

6.5 



4 

4 


BN12 

PTFlOufj 

-9.0 

4.1 



4 

4 


BN12 

PTFlOwnm 

-7.1 

4.8 



4 

4 


BN12 

PTFlOwof 

-3.9 

9.1 



4 

4 


BN12 

PTFlOxyt 

-4.6 

7.3 



2 

5 


BN12 

SN2011fe 

-16.0 

- • . 45.9 



32 

35 

32 

M12 

SN2014J 

-10.0 

- • • 72.4 



24 

24 

24 

F14 


References: MOO - Meikle (2000); Ph06 - Phillips et al. 

(2006); K04a - Krisciunas et al. (2004); VOS - Valentini et 
al. (2003); K04b - Krisciunas et al. (2004); P08 - Pignata et 
al. (2008); CIS - Cartier et al. (2014); ER06 - Elias-Rosa et 
al. (2006); K09 - Krisciunas et al. (2009); P07 - Pastorello 
et al. (2007); CSP - Carnegie Supernova Project Contreras 
et al. (2010); Stritzinger et al. (2011); M12 - Matheson et 

al. (2012); F14 - Foley et al. (2014) 
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We have included in our sample near-IR SN la photom¬ 
etry from Meikle (2000) and several SNe la observed by the 
European SN Consortium (ESC; Benetti et al. 2004; Pig- 
nata et al. 2008; Elias-Rosa et al. 2006; Pastorello et al. 
2007; Krisciunas et al. 2009). Twelve SNe la have been dis¬ 
cussed by Barone-Nugent et al. (2012) and have data only 
near the first maximum. We also included NIR photometry 
from two recent nearby explosions, SN2011fe (Matheson et 
al. 2012) and SN2014J (Foley et al. 2014). 

The 91 objects used in this work are listed in Table 1, 
where the phase range of observations (first and last ob¬ 
servation), total number of observations in each filter and 
the reference for each data set are tabulated. The sample is 
dominated by SNe la from the CSP and we show the results 
separately for the CSP and non-CSP objects. It is worth 
noting that there are 15 SNe la with observed NIR light 
curves beyond 100 days. 

As can be seen in Table 1 and displayed in Fig. 1, the 
A-band light curves are sparsely sampled and not enough 
objects are available for detailed analysis. Therefore, we ex¬ 
clude the A-band light curves from further analysis. 


3 NIR LIGHT CURVE MORPHOLOGY 

In parameterising the NIR light curves we follow the nomen¬ 
clature introduced by Biscardi et al. (2012). The first maxi¬ 
mum in filter X, Mi{X) is reached at a phase, ti{X) relative 
to the phase of the B maximum d). The light curves 

dip to a minimum Mq{X) at to{X) before reaching a second 
maximum, M 2 {X) at time, t 2 {X). These three prominent 
features are shown in Fig. 1. The light curves are plotted 
without normalisation for phase or corrections for possible 
differences in photometric systems or absorption. 

3.1 Light curve fitting 

We fit the optical light curves using the programme SNooPy 
(Burns et al. 2011) to determine the peak of the R-band 
light curve through cubic-spline fitting and determine the 
phase, used throughout this work and the peak bright¬ 
ness . SNooPY also determines the Amis parameter^ 

commonly used to characterise the SN light curve shape and 
provides an estimate of the extinction in the host galaxy. The 
SNooPy light curve fit parameters and distance moduli for 
the SNe la in our sample are given in Table 2. We used the 
published values of the distance modulus, fj,, for non-CSP ob¬ 
jects, (references in Table 2). The distance moduli for CSP 
SNe la not in the Hubble flow are taken from Contreras et 
al. (2010) and Stritzinger et al. (2011), and the individual 
references are listed in Table 2. The distance moduli for the 
rest of the SNe are based on the host galaxy redshift from 
the NASA/IPAC Extragalactic Database adopting a Hubble 
constant of Hq=7Q km s“^ Mpc“^. 

While the extinction is much reduced in the NIR it 
is not entirely negligible. Using the Cardelli extinction law 
(Cardelli, Clayton & Mathis 1989), the extinction in the H- 
band is ~18% of that in V-band. We have also included 


Table 2. Timing of Bmax, Amis, distance modulus 
for the SN in our sample. References for the published 
distance moduli are provided in the footnote. 


SN Name 

tiBmax) 

(MJD) 

err 

Amis 

(mag) 

err 

M 

(mag) 

err 

SN1980N 

44585.8 

0.5 

1.28 

0.04 

31.59 

0.10 

SN1981B 

44672.0 

0.2 

1,10 

0.04 

30,95 

0,07 

SN1986G 

46561.0 

0.1 

1.76 

0.10 

28.01 

0.12 

SN1998bu 

50953.3 

0.5 

1,01 

0.02 

30,20 

0,10 

SN1999ac 

51251.0 

0.5 

1.34 

0.02 

33.50 

0.30 

SN1999ee 

51469.1 

0.5 

1.09 

0.02 

33.20 

0.20 

SN1999ek 

51481.8 

0.1 

1.17 

0.03 

34.40 

0,30 

SN2000E 

51557.0 

0.5 

0.99 

0.02 

31,90 

0.40 

SN2000bh 

51635.2 

0.5 

1.16 

0.01 

34,60 

0,30 

SN2001ba 

52034.5 

0.5 

0.97 

0.05 

35.40 

0.50 

SN2001bt 

52063.4 

0.5 

1.18 

0.02 

34.10 

0,40 

SN2001cn 

52071.6 

0.2 

1.15 

0.02 

34.10 

0.30 

SN2001CZ 

52103.9 

0.1 

1.05 

0.07 

33.50 

0,10 

SN2001el 

52182.5 

0.5 

1.13 

0.04 

31.30 

0.20 

SN2002bo 

52356.5 

0.2 

1.12 

0.02 

31.80 

0,20 

SN2002dj 

52450.0 

0.7 

1.08 

0.02 

32.90 

0.30 

SN2002fk 

52547.9 

0.3 

1.02 

0.04 

32.59 

0.15 

SN2003cg 

52729.4 

0.5 

1.12 

0.04 

31.28 

0,20 

SN2003hv 

52891.2 

0.3 

1.09 

0.02 

31,40 

0.30 

SN2004ef 

53264.4 

0.1 

1.45 

0.01 

35.57 

0,07 

SN2004eo 

53278.4 

0.1 

1.32 

0.01 

34.03 

0.10 

SN2004ey 

53304.3 

0.1 

1.02 

0.01 

34,01 

0,12 

SN2004gs 

53356.2 

0.1 

1.53 

0.01 

35.40 

0.08 

SN2004gu 

53362.2 

0.2 

0.80 

0.01 

36.59 

0,04 

SN2005A 

53379.7 

0.2 

1.08 

0.02 

34.51 

0.11 

SN2005M 

53405.4 

0.1 

0.80 

0.04 

35.01 

0,09 

SN2005ag 

53413.7 

0.2 

0.87 

0.01 

37.80 

0,03 

SN2005al 

53430.5 

0.1 

1.30 

0.01 

33.79 

0.15 

SN2005am 

53436.9 

0.1 

1.48 

0.01 

32,85 

0,20 

SN2005el 

53647.0 

0.1 

1.40 

0.01 

34.04 

0.14 

SN2005eq 

53654.4 

0.1 

0.82 

0.01 

35.46 

0,07 

SN2005hc 

53666.7 

0.1 

0.80 

0.01 

36.50 

0.05 

SN2005hj 

53673.8 

0.2 

0.80 

0.02 

37,03 

0,04 

SN2005iq 

53687.7 

0.1 

1.28 

0.01 

35.80 

0.15 

SN2005kc 

53697.7 

0.1 

1.12 

0.02 

33.89 

0,15 

SN2005ki 

53705.5 

0.1 

1.36 

0.01 

34.73 

0.10 

SN2005na 

53740.2 

0.1 

1.03 

0.01 

35.34 

0,08 

SN2006D 

53757.7 

0.1 

1.47 

0.01 

33.00 

0,15 

SN2006X 

53786.3 

0.1 

1.09 

0.03 

30.91 

0.08 

SN2006ax 

53827.2 

0.1 

1.04 

0.01 

34.46 

0,11 

SN2006bh 

53833.6 

0.1 

1.42 

0.01 

33.28 

0.20 

SN2006br 

53853.7 

0.6 

1.45 

0.05 

35.23 

0,08 

SN2006ej 

53976.4 

0.2 

1.37 

0.01 

34.62 

0.11 

SN2006eq 

53975.9 

0.4 

1.88 

0.04 

36.66 

0,04 

SN2006et 

53993.7 

0.1 

0.88 

0.01 

34.82 

0.10 

SN2006ev 

53990.1 

0.3 

1.34 

0.01 

35.40 

0,08 

SN2006gj 

54000.3 

0.2 

1.56 

0.04 

35.42 

0.08 

SN2006gt 

54003.1 

0.3 

1.71 

0.03 

36.43 

0,05 

SN2006hb 

54006.0 

0.3 

1.69 

0.02 

34.11 

0,13 

SN2006hx 

54021.9 

0.2 

1.07 

0.05 

36.47 

0.05 

SN2006is 

54007.5 

0.4 

0.80 

0.01 

35.69 

0,07 

SN2006kf 

54041.3 

0.1 

1.51 

0.01 

34.78 

0.10 

SN20061U 

54034.4 

0.2 

0.92 

0.01 

36,92 

0,04 

SN2006ob 

54063.4 

0.1 

1.51 

0.01 

37.08 

0.04 

SN2006OS 

54063.9 

0.2 

1.08 

0.02 

35.70 

0,07 

SN2007A 

54113.1 

0.2 

1.06 

0.04 

34,26 

0.13 

SN2007S 

54143.8 

0.1 

0.81 

0.01 

34.06 

0,14 

SN2007af 

54174.4 

0.1 

1.11 

0.01 

32.10 

0.10 

SN2007ai 

54173.5 

0.3 

0.84 

0.02 

35.73 

0.07 

SN2007as 

54181.3 

0.4 

1.27 

0.03 

34.45 

0,12 

SN2007bc 

54200.3 

0.2 

1.27 

0.02 

34.89 

0.10 

SN2007bd 

54206.9 

0.1 

1.27 

0.01 

35.73 

0,07 

SN2007bm 

54224.1 

0.2 

1.11 

0.02 

32.30 

0.07 

SN2007ca 

54227.7 

0.2 

1.05 

0.03 

34,04 

0,14 

SN2007if 

54343.1 

0.6 

1.07 

0.03 

37.59 

0.03 

SN2007jg 

54366.1 

0.3 

1.09 

0.04 

36.03 

0,06 

SN20071e 

54399.3 

0.1 

1.03 

0.02 

32.34 

0.08 

SN2007nq 

54398.8 

0.1 

1.49 

0.01 

36.44 

0,05 

SN2007on 

54419.8 

0.4 

1.65 

0.04 

31.45 

0.08 

SN2008C 

54466.1 

0.2 

1.08 

0.02 

34.34 

0.12 

SN2008R 

54494.5 

0.1 

1.77 

0.04 

33.73 

0,16 

SN2008bc 

54550.0 

0.1 

1.04 

0.02 

34.16 

0.13 

SN2008bq 

54562.1 

0.2 

0.78 

0.02 

35.79 

0,06 

SN2008fp 

54730.9 

0.1 

1.05 

0.01 

31.79 

0.05 

SN2008gp 

54779.1 

0.1 

1.01 

0.01 

35.79 

0,06 

SN2008hv 

54817.1 

0.1 

1.30 

0.01 

33.84 

0.15 

SN2008ia 

54813.2 

0.1 

1.34 

0.01 

34,96 

0,09 

SN2011fe 

55815.0 

0.3 

1.20 

0.02 

28.91 

0.20 

SN2014J 

56689.7 

0.3 

1.10 

0.02 

27.64 

0,10 


References: SN1980N, SN1981B, SN1986G, 

SN1998bu (Meikle 2000); SN1999ac, SN1999ee, 

SN2002dj (Willick et al. 1997); SN2000E, SN2002bo 
(Tully 1988); SN2001el (Ajhar et al. 2001); SN2003hv, 
SN2007on (Tonry et al. 2001); 2006X (Freedman et 
al. 2001) (note: we add the correction to p’s from 
Tonry et al. (2001) and Ajhar et al. (2001) using the 
values provided in Jensen et al. (2003)) 


^ The Amis calculated by SNooPy from all available filters. It 
is linearly related to Amis(B) (see Burns et al. 2011) 
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Phase{d) Phase(d) 


Figure 1. The Y, J, H and K absolute magnitude light curves of SNela. The CSP sample is shown as blue pentagons. The green 
triangles indicate the non-CSP sample (see Table 1 for the references). The root-mean-square (RMS) scatter for each day is shown in 
the lower panels. For guidance a line at a scatter of 0.2 magnitudes is drawn in each panel. The luminosities are based on the literature 
distances (cf. Table 2) and no correction has been applied. The steady increase of the scatter after the first maximum is evident. 


some heavily-extinguished SNela like SNe 1986G, 2005A, 
2006X, 2006br and 2014J without an extinction correction 
in our sample. Therefore, the observed scatter is larger than 
the intrinsic variation among SNela. We have chosen not to 
apply a correction for the host galaxy extinction as the red¬ 
dening law remains under debate (e.g. Phillips et al. 2013). 

We fit a spline interpolation to the data to derive the 
phase and magnitude at maximum, the minimum and the 
second maximum in each filter. In order for a measurement 
of the minimum and second maximum to be made, we re¬ 
quire ^4 observations at late phases (>7 d for the minimum 
and >15 d for the second maximum). We also require that 
observations at least 4 d before tmax were available. The 
uncertainties for each derived parameter were calculated by 


repeating the fits to 1000 Monte Carlo realisations of light 
curves generated using the errors on the photometry. 

The NIR light curves are very uniform up to the time of 
the R-band maximum, ~3-4 d after the maximum is reached 
in the NIR light curves. In the lower panels in Fig. 1, we show 
the RMS scatter for each epoch. 

The scatter remains small for ~1 week around maxi¬ 
mum in the different bands independent of the sample (CSP 
or non-CSP). The two samples show very consistent scatter 
out to a phase of ~10 d in J and ~25 d in H, after which 
they start to deviate. With only 14 objects in the non-CSP 
sample in this phase range compared to 25 from the CSP, 
we consider the differences not statistically significant. The 
scatter continues to increase beyond 35 d, which we attribute 
to the colour evolution (see §3.6 below). 
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Figure 2. Distribution of ti for the Y, J, and H light curves 
relative to the S-band maximum. The filled histograms are for 
SNe la from the CSP sample while the open histograms show the 
combined CSP and non-CSP sample. In the V band, we only use 
the objects observed by the CSP. The NIR light curves peak at 
within —2 to —7 d relative to the B maximum. 


3.2 The first maximum 

Elias et al. (1981) showed, in a small sample of SNe la, that 
the JHK light curves of SNe la peak earlier than the optical 
light curves. We confirm this for our sample (Fig. 2). 

The NIR light curves peak within —2 to —7 d of the B- 
band peak confirming the result of Folatelli et al. (2010) for 
SNe la with Amis (i3)< 1.8. There is no obvious difference 
between our full sample and the CSP sample as seen in the 
scatter. The distribution of ti is remarkably tight (less than 
one day dispersion) for all filters. This indicates a close rela¬ 
tionship between and the NIR ti values for the SNe la 
in our sample. 

Table 3 gives the phase of lowest scatter measured in 
each filter. Without any attempt to normalise the light 
curves, we find the smallest scatter in all NIR light curves 
near ti. The dispersion remains very low for ~1 week, before 
increasing to >0.2 mag at later phases. 


3.3 The minimum 

The minimum in J occurs ~2 weeks after (Fig- 3). 

The Y light curves dip about 3 d earlier at to = -1-11 d. 


Figure 3. Distribution of to in Y, J, and H light curves. The 
Y and H minimum is reached a few days before the J light curve 
dips, which occurs ~2 weeks after B-band maximum. 


Table 3. Magnitude scatter in NIR light curves. The second col¬ 
umn indicates the time of minimum magnitude scatter in each 
filter, while the third column gives the scatter at this phase. The 
fourth column gives the phase range for which the scatter stays 
below 0.2 mag. 


Filter 

t 

(days) 

cf{M) 

(mag) 

Phase range 
(f7(M) < 0.2 mag) 

SN sample 

Y 

-4.4 

0.15 

[-4 , +1] 

CSP 

J 

-3.6 

0.16 

[-4 , +3] 

CSP 

J 

-3.8 

0.17 

[-6 , +1] 

non-CSP 

H 

-5.1 

0.17 

[-5 , +1] 

CSP 

H 

-4.7 

0.14 

[-7 , +2] 

non-CSP 


The minimum in H is reached on average about 2 d before 
J at to = -1-12 d. The phase range is still relatively narrow 
with the minima all occurring within roughly ±2 d. While Y 
and H display a tight distribution of to, the J distribution 
exhibits a tail of late minima. 

No significant difference between the CSP and the liter¬ 
ature samples can be seen in the distributions. The scatter 
among Mo is fairly small for the three NIR bands although 
about 2 to 3 times larger than at ti. In the y-band, the 
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Figure 4. Distribution of t 2 in the Y, J and H light curves. 
Note the expanded scale on the abscissa compared to the phase 
ranges in Figs. 2 and 3. 

scatter is low (<0.2 mag) immediately after to, while the J 
and H light curves display a larger dispersion at this point. 

3.4 The second maximum 

Figure 4 shows that to occurs over a wide range of phases 
and can vary by as much as 20 d from one SN la to another. 
This diversity had been observed before for the i light curves 
(e.g. Hamuy et al. 1996; Folatelli et al. 2010) and for JHK 
(Mandel et al. 2009; Biscardi et al. 2012). In Figure 4, we 
show that the mean to is later in T- and J-bands, with 
the Tf-band light curves reaching to, on average, a few days 
earlier. The scatter in the luminosity starts to increase slowly 
after to, and is seen to increase significantly (>0.5 mag) 
around the time of to (Fig. 1). 

3.5 The late decline 

After the second maximum, the light curves steadily decline. 
In Fig. 5 we show the distribution of slopes calculated for 
SNe la with at least three observations in phases 40<t<90 
d. This choice of phase range ensures that the measurements 
are not influenced by the second maximum. All the SNe la 
in our sample with data at these late phases come from the 
CSP. 


Figure 5. Distribution of decline rates after to- 

There is a tight distribution of decline rates in the Y- 
and 7f-bands, with the notable exception of SN 2005M in H, 
which declined twice as fast as the other SNe la. The scatter 
in the 77-band decline rate, when excluding SN 2005M, is 
only 0.004 mag per day, similar to the scatter in the Y filter. 
The same SN is also one of the two fast declining objects in 
J (in the fastest bin, middle panel), where the scatter in the 
decline rates is significantly larger than for the other two 
bands. 


3.6 NIR colours 

Elias et al. (1985) showed that the early J — 77 colour evo¬ 
lution is rather uniform for SNe la. In Fig. 6, we show the 
colour curves (Y — J and J — 77) for the CSP and non-CSP 
objects. The scatter at each epoch is plotted in the lower 
panels. Similar to the light curves, the early colour evolu¬ 
tion is similar for most SNe la in our sample. 

At the first maximum the scatter is minimal inY — J 
{a(Y — J) = 0.07 mag) and inj —77 {a{J — H) — 0.05 mag). 
The scatter stays <0.1 mag between —7 and —3 d for Y — J 
and between —10 and —3 d for J — 77. At early times, we 
find that the samples display the same scatter. Immediately 
after the first maximum, the colour curves start to deviate 
like the light curves. By the time of the minimum the colours 
display a wide spread which continues to increase into the 
late decline. 
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Figure 6. Left-. Y — J colour curve. The mean ti, to and to in the Y filter are plotted as green lines. The scatter plot in the lower panel 
only contains information from the CSP objects since the non-CSP sample does not have F-filter observations. The blue points are the 
objects in the CSP sample, whereas the green points are the non-CSP objects. Right: J — H colour curve. The epoch of minimum scatter 
in J, the average value of to and to in J are overplotted in green. The lower panel shows the evolution of the scatter around the mean 
in the colour curve. 


The Y — J colour remains fairly constant until about one 
week after maximum when the Y — J colour starts a mono¬ 
tonic evolution towards bluer colours. For the J — H colour, 
SNela evolve to redder colours after maximum. Then after 
the light curve minimum, the J — H colour tends towards 
slightly bluer colours until ~30 d. 


4 CORRELATIONS 

The uniformity of the NIR light curves at maximum light 
suggests that the size of the surface of last scattering at these 
wavelengths is independent of the details of the explosion 
and progenitor. 

A diversity in the NIR only becomes visible at later 
phases as the lines contributing to the line blanketing opac¬ 
ity arise from deeper in the ejecta. The phase and magnitude 
the minimum and in particular the second maximum of the 
NIR light curves display large variations. Correlating these 
changes with other SN parameters should shed light on the 
physical processes underlying the explosions and the release 
of the radiation. 

In Fig. 7 we show example NIR light curves extending 
to late times for 10 SNe la in our sample. It should be noted 
that among these 10 SNe are two where the extinction has 
been determined to be high (SN 2006X and SN 2008fp). 
Therefore, some of the scatter can be attributed to uncor¬ 
rected host galaxy extinction. 

In the following sections, we note that correlations re¬ 


ported with r > 0.4 are significant, and those with r > 0.65 
are termed as strong. 

4.1 NIR light curve properties 

In the previous section, we found a large diversity in the 
NIR light curve properties at post-maximum epochs. In this 
section, we investigate correlations between these properties 
in more detail. We find a correlation between to and to in JH 
filters (Pearson parameter r ~ 0.5) suggesting that a later 
minimum is followed by a later second maximum. We also 
find a correlation between Mo and Mo in Y and J filters, 
with r = 0.59 and r = 0.50 respectively. This implies that 
SNe la with a more luminous minimum also display a more 
luminous second maximum. 

Interestingly, the Mo and to do not correlate (r < 0.4) in 
any filter. However, there we do find a significant correlation 
across different filters between Mo and to (e.g. Mo in Y band 
correlates with to in J band with r = 0.61). Mandel et al. 
(2009) found that the rate of luminosity increase (rise rate 
rjP in the Mandel et al. 2009 nomenclature) to the second 
maximum correlated with the luminosity of the first peak. 
We cannot confirm this at any significance with our data. 

Comparison of Figs. 3 and 4 reveals that the light curve 
evolution in the T-band is slowest amongst the NIR filters. 
SNe la on average reach to in Y earlier than in J and H, 
but reach to later in the Y compared to the other filters. 
The rise time in Y is nearly 4 d longer than in J and nearly 
7 d longer than in H. 

The phase and luminosity of the second maximum 
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Figure 7. Complete absolute magnitude light curves in YJH 
for the 10 objects for which a late decline is measured in all three 
filters. The figure illustrates that the decline rates at late times 
are uniform, whereas there is a large scatter in the absolute mag¬ 
nitude at +55 d indicated by the vertical line. 


strongly correlate between the NIR filters (Fig. 8 ). SNe la 
with a later t 2 display a higher luminosity during the late 
decline. The luminosity at 55 d past the B maximum, here¬ 
after referred to as Mjss was chosen to ensure all SNe la 
have entered the late decline well past t 2 . A choice of a later 
phase may be more representative of the decline but would 
result in a smaller SN sample as not many objects are ob¬ 
served at these epochs and the decreasing flux results in 
larger uncertainties. 

In Figure 9 we plot Mjss against t 2 and M 2 . A clear 
trend between MI 55 and t 2 is present in all filters (Pearson 
coefficients r of 0.78, 0.92, 0.68 in the YJH, respectively; 
Fig. 9, left panels). At this phase SNe la are on average 
most luminous in Y followed by H and J, a trend that is 
already present at M 2 . This is also reflected in the NIR 
colour evolution (see section §3.6). 

The dispersion in MI 55 is large with (j(M| 55 ) = 

0.48 mag, 0.51 mag and 0.30 mag in Y, J and H, respec¬ 
tively. This is not unexpected as it continues the trend to 
larger (luminosity) differences in the NIR light curves with 
increasing phase. 



Figure 8. Phases (left) and luminosity {right) of the second 
maximum in NIR filters. There are clear correlations between the 
filters with the weakest trend in the H versus J luminosities. The 
black line is a one-to-one relation. 



10 20 30 -17 -18 -19 

t^[d) M^imag) 


Figure 9. (Left:) Absolute magnitude at t = 55 d in Y JH vs. t 2 . 
(Right:) MI 55 compared to the absolute magnitude of the second 
maximum M 2 - 
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Figure 10. Comparison of t 2 in the NIR light curves with 
Ami5(S). 


Figure 11. M 2 versus Amis is shown and a weak correlation is 
observed. 


4.2 Correlations with optical light curve shape 
parameters 

It is interesting to see whether the NIR light curve parame¬ 
ters correlate with some of the well known optical light curve 
shape features (Amis; Burns et al. 2011). 

Folatelli et al. (2010) have shown that the value of t 2 
in i correlated with Amis(i3). Since the dispersion increases 
with phase we concentrate on the second maximum and ex¬ 
plore whether the timing and the strength correlate with the 
optical light curve shape. 

The second maximum in the NIR is a result of an 
ionization transition of the Fe-group elements from dou¬ 
bly to singly ionized states. Models predict that t 2 depends 
on the amount of Ni (Mni) synthesized in the explosion 
(Kasen 2006). Since Mjvi is known to correlate with 
which itself is a function of the light curve shape (Arnett 
1982; Stritzinger et al. 2006; Mazzali et al. 2007; Scalzo et 
al. 2014), we explored the relation between t 2 and Amis 
(Tab. 2). 

Figure 10 confirms a strong correlation between t 2 and 
Amis. The Pearson coefficients are 0.91, 0.93 and 0.75 for 
the Y, J, and JT-bands, respectively. The weaker correlation 
in H may be due to a low dependency on t 2 for objects with 
Amis < 1-2, where t 2 appears independent of Amis- This 
applies only to SNe la with a slow optical decline. For higher 
Amis, the trend is as strong as in the other filters. 


A linear regression yields equation 1 for each filter. The 
RMS scatter is 2.1 d, 1.8 d and 3.0 d in the YJH filters. 

t 2 {Y) = (-20.6 ± 1.0)Amis -b (53.6 ± 1.2) (la) 

t 2 {J) = (-20.3 ± 0.9)Amis -b (51.7 ± 1.1) (lb) 

t 2 {H) = (-16.0 ± 2.7)Amis -b (42.3 ± 2.8) (Ic) 

Since the phase of the second maximum ( 12 ) in NIR 
bands is strongly correlated with Amis and thereby the op¬ 
tical maximum luminosity, t 2 could also serve as an indicator 
of the luminosity of SNe la. It is noteworthy that even ex¬ 
treme cases, like SN 2007if, which have been associated with 
super-Chandrasekhar mass progenitors (Scalzo et al. 2010, 
2012) are fully consistent with the derived relations. 

There is a weak correlation between M 2 and Amis as 
shown in Fig. 11, with very little difference between objects. 
We find r values of 0.59, 0.50 and 0.63 for the Y, J, and H 
filters, respectively. 


5 DISCUSSION 

As has been shown earlier in this paper and in a num¬ 
ber of other publications (Meikle 2000; Wood-Vasey et al. 
2008; Barone-Nugent et al. 2012; Weyant et al. 2014), SNIa 
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NIR light curves show remarkable uniformity around max¬ 
imum light, compared to at optical wavelengths. This uni¬ 
formity, combined with reduced effects of extinction, holds 
great promise for the use of SNe la as distance indicators in 
the NIR. Kattner et al. (2012) proposed to further reduce 
the scatter in the luminosity of the first maximum in the 
NIR by a decline rate correction similar to the procedure in 
the optical. 

The spread in optical luminosity is attributed to differ¬ 
ent masses, or distributions, of ®®Ni within the ejecta (Ar¬ 
nett 1982; Stritzinger et al. 2006; Scalzo et al. 2014). The 
indifference of the NIR maximum light to the nickel mass 
suggests that it is intermediate-mass elements that domi¬ 
nate the opacity in these bands at maximum light. 

The NIR light curves showed an increased dispersion at 
later times. We attribute this increased scatter to differences 
in the speed of the evolution of the SNe la. The phase of the 
second maximum depends on the mass and distribution of 
®®Ni, the change in opacity, the ionisation and the dominant 
species setting the emission. 

Not all SNe la display a second maximum (e.g. Krisciu- 
nas et al. 2009) and we restricted our analysis only to objects 
in which the second maximum is dehned well enough to be 
ht. This translates into a sample including only SNe la with 
Ami 5 (R) < 1.8 mag. Events without a second maximum 
tend to be of low luminosity, often similar to SN 1991bg, 
and with large 

The strength of the second peak in the iJHK light 
curves does not correlate with Amis, but the phase does 
(Folatelli et al. 2010; Biscardi et al. 2012). 

As presented in the previous section, we conhrm this 
relation between t 2 and Amis for the Y JH filters. The cor¬ 
relation of M 2 with Amis is rather weak (Fig. 11) although 
it appears somewhat stronger than in Biscardi et al. (2012) 
(the Pearson coefficients in our data are 0.5 and 0.63 for J 
and H compared to 0.12 and 0.08 in Biscardi et al. 2012). 

5.1 A possible physical picture 

The various features of the NIR light curves can be assem¬ 
bled into a physical picture of the explosions. The striking 
similarities of the late decline rates, when the SN becomes 
increasingly transparent to the 7 —rays generated by the ra¬ 
dioactive decays, indicate that the internal structure of the 
explosions is probably similar for all SNe la considered here. 
The uniform decline rates are consistent with the predictions 
of Woosley et al. (2007) for a range of Chandrasekhar-mass 
models, with different Mjvi but similar radial distribution of 
iron-group elements. We find that the late-time decline rate 
in J is faster than in Y and H, a trend also seen between the 
simulated J and H light curves in Dessart et al. (2014). The 
NIR light curves depend very little on the explosion geom¬ 
etry (Kromer & Sim 2009). The DDC15 models of Blondin 
et al. (2013) show that the YHK decline rates are similar 
to the pseudo-bolometric decline rate (S. Blondin, private 
communication). The J band shows a faster decline due to 
a lack of emission features (Spyromilio, Pinto, & Eastman 
1994). This also explains the evolution of the J — H colour 
curve to redder colours at late times. 

At the these late times MI 55 shows a large scatter (cf. 
Fig. 9). If the similar late decline rates ( cf. Fig. 5) indicate 
a similarity in the evolution of the 7 -ray escape fraction for 





Figure 12. The phase of the second maximum versus tj^. The 
black line is the one-to-one relation. 


different SNe la,then, a higher luminosity would translate 
into a higher energy input at late phases, i.e. a larger mass 
of Fe-group elements. 

Kasen (2006) predicted that the second maximum 
should be delayed for larger Fe masses, which is exactly 
what is observed in the NIR (see also Jack, Hauschildt, & 
Baron 2012). According to Kasen & Woosley (2007), the 
faster decline in the B-band light curve is mostly due to 
line blanketing through Fe and Co lines, which shifts the 
emission into the NIR and shapes the NIR light curves after 
maximum. The optical colour evolution post B-maximum 
is suggested to be more rapid for explosions with lower Ni 
masses. If this is true then the onset of the uniform B — V 
colour evolution (referred to as the ’Lira law’ and originally 
defined as the uniformity of the slope of the colour curve 
from 30 to 90 days, although the onset is generally at ear- 
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Figure 13. A summary of the correlations between the timing parameters (ti, to and to), Amis i iL and the late time luminosity 
(Miss), ia each filter There is a strong correlation amongst the t 2 , Amis, and M|ss (r > 0.65). The colorbar scales from white to 
red in ascending order of correlation strength. Correlations with r < 0.4 have been set to white 


lier phase; Phillips et al. 1999) marks the beginning of the 
nebular phase. At these epochs, the SN la spectrum is dom¬ 
inated by emission lines from Fe-group elements and the 
emission line strength depends on the Fe mass in the explo¬ 
sion. We measure the time at which the SN enters the Lira 
law, hereafter tn, as the epoch of inflection in the (B — V) 
colour curve, at late times. The procedure for measuring tr 
is identical to the measurement of t 2 described in §3.4. 

Figure 12 shows that t 2 coincides nearly exactly with 
tL for the Y and J bands. While the H light curve peaks 
3-4 d earlier. The reason for this is not entirely clear, but 
both Y and J-bands are expected to be strongly influenced 
by Fe lines, while H is dominated by Co lines (e.g. Marion 
et al. 2009; Jack, Hauschildt, & Baron 2012). The striking 
coincidence of t 2 and II in the NIR light curves is further 
evidence that t 2 directly depends on the Ni mass in the 
explosion. 

Scalzo et al. (2014) found that the Ni mass depends on 
Ami 5 (B) and the bolometric peak luminosity. All these pa¬ 
rameters , Amio{B), t 2 and tL, and NIR late-phase 

luminosity) can be tied to the Ni (and/or Fe) mass in the 
explosion. We provide a summary of the important corre¬ 
lations analysed in this work for the three NIR filters in 
Figure 13, where we correlate the timing parameters {ti, to, 
t 2 ), Amis, tL and NIR late phase luminosity (Mjss). With 
most of the Fe synthesised in the explosion, we propose that 
all light curve parameters point to the Ni as the dominant 
factor in shaping the light curves. We find a consistent pic¬ 
ture that the properties of second maximum in the NIR light 


curves are strongly influenced by the amount of Ni produced 
in the explosion followed by a more luminous decline. 

The NIR colours show a pronounced evolution after ti, 
the flux in the J band decreases significantly with respect 
to both Y and H shortly after ti (Fig. 6). The reason for 
this depression is most likely the lack of transitions provid¬ 
ing the required channel for radiation to emerge around 1.2 
/rm (Spyromilio, Pinto, & Eastman 1994; Hoflich, Khokhlov 
& Wheeler 1995; Wheeler et al. 1998). This persists until 
t 2 when the [Fe II] (1.25 pm) emission line forms and the 
J-band magnitude recovers relative to H, although not com¬ 
pared with Y, which is also dominated by Fe lines (Marion 
et al. 2009). The Bt filter is dominated by Co II lines (e.g. 
Kasen 2006; Jack, Hauschildt, & Baron 2012). After the sec¬ 
ond maximum J — H turns redder again due to the faster 
decline rate in the J filter. 

We investigated whether the (optically) fast-declining 
SNela can be grouped within our analysis. These SNe have 
very low ®®Ni mass which is almost certainly centrally lo¬ 
cated (low velocities of Fe III in the optical at late times in 
SN 1991bg; Mazzali et al. 1998) and transition so rapidly 
that a second maximum barely has time to form. Thus, we 
could not include them directly in our study. However, a 
comparison of M\ss with Amis{B) can be made for four ob¬ 
jects of this class (SNe 2005ke, 2006mr, 2007N and 2007ax). 
The trend of these objects showing a fainter Mjss luminosity 
with larger Amis{B) is followed, but it is unclear whether 
these SN 1991bg-like SNela follow the same relation as their 
brighter counterparts. The result remains inconclusive as the 
scatter remains currently too large. We note that the decline 
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rates in the NIR at late times for these objects do not differ 
signihcantly from our sample. They follow the same distri¬ 
bution as in Fig 5. 

5.2 Improved Distance Measurements? 

The uniformity of the NIR light curves is in use for cos¬ 
mological projects (Barone-Nugent et al. 2012; Weyant et 
al. 2014). Krisciunas et al. (2009) hnds no correlation of the 
maximum luminosity with other light curve parameters, but 
identihed the (optically) fast-declining SNela as sublumi- 
nous in the NIR compared to the other SNela in their sam¬ 
ple. In Kattner et al. (2012), the authors hnd a weak trend 
between the NIR first maximum luminosity and Amis. 

We looked for NIR light curve properties to further im¬ 
prove SNela as distance indicators. In Fig. 14 we find no 
signihcant evidence for a correlation between Mi and t 2 for 
objects in our sample. The faintness of SN 2006X in this 
diagram is probably due to the strong absorption towards 
this SN. 

The scatter in Fig. 1 falls to around 0.2 mag at later 
NIR phases. Between to and t 2 all NIR light curves reach 
comparable luminosities. These phases are between 10 and 
20 d in Y, near 15 and 20 d in J and H when SNe la are only 
about 1 (0.5) magnitude fainter than at the first maximum 
in Y (H). The decline in the J light curves is quite steep after 
the hrst maximum and the SNe have already faded by nearly 
2 magnitudes. At least for Y and H, it might be worthwhile 
to investigate whether good distances can be determined at 
later phases (between +10 and +15 days). The advantage 
would be a distance measurement with a reduced extinction 
component, i.e. mostly independent of the exact reddening 
law, and the possibility of targeted NIR observations even 
when the first maximum had been missed. 


6 CONCLUSIONS 

The cosmological interest in SNe la in the NIR stems from 
the small observed scatter in the peak magnitudes. We con- 
Hrm this with our extended literature sample despite our 
simple assumptions on distances and neglect of host galaxy 
absorption. The phase of the hrst maximum in the YJH 
light curves shows a narrow distribution. The uniformity of 
the SNe la in the NIR lasts until about one week past B 
maximum. The NIR light curves diverge showing a large 
scatter by the time of the second maximum and thereafter. 
The IR colour curves are uniform at early phases with in¬ 
creasing scatter after 20 days. 

These hndings corroborate the use of only few IR obser¬ 
vations near the hrst maximum to obtain good cosmological 
distances (Barone-Nugent et al. 2012; Weyant et al. 2014). 
The small scatter prevails for nearly a week around the hrst 
maximum and potentially another window opens near the 
second maximum, where the scatter again appears rather 
small. A condition for these measurements is accurate phase 
information to be able to use sparse NIR observations to de¬ 
rive the distances. The phase information could come from 
accurate optical light curves, as used here, or through spec¬ 
tra. The latter are needed for redshifts and classihcation in 
any case. 

The information contained in the NIR light curves 



Figure 14. The magnitude of the first maximum in the Y JH 
bands is plotted against the phase of the second maximum in J 
band. The underluminous object in Y and J bands is the heav¬ 
ily extinguished SN 2006X. In the J band SN 2014J also appears 
fainter. 

points towards the nickel/iron mass as the reason for the 
variations. The absolute magnitude 55 d past maximum to¬ 
gether with the very uniform decline rate of the light curves 
in all bands, including the optical (Barbon, Ciatti, & Rosino 
1973; Phillips et al. 1999; Leibundgut 2000), points towards 
differences of the energy input into the ejecta with a rather 
uniform density structure. This is also found in theoreti¬ 
cal studies (Kasen & Woosley 2007; Kromer & Sim 2009; 
Dessart et al. 2014), although most of these models employ 
Chandrasekhar-mass explosions and a wider range of ejecta 
mass models may need to be explored to confirm the simi¬ 
larities in the structure of the ejecta. A higher luminosity at 
a fixed phase in the nebular phase points towards a larger 
iron core and hence a higher nickel production in the explo¬ 
sion. A corollary is the phase of the second maximum which 
also occurs later for larger iron cores (Kasen 2006). 

The strong correlation of the phase of the second IR 
maximum with the optical light curve shape parameter 
Amis and the onset of the uniform [B — V) colour evolu¬ 
tion (‘Lira law’) point towards SNe la as an ordered family 
and nickel mass as the dominating factor in shaping the ap¬ 
pearance of SNe la. With a higher nickel mass a larger Fe 
core is expected. This would result in higher expansion ve¬ 
locities observed at late phases (Mazzali et al. 1998). It is 
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worth checking whether this prediction holds true in future 
observations. 

The phase of the second maximum should provide a 
handle to determining the nickel masses in SN la explosions, 
in particular for SNe la where absorption is significant. In 
cases like SN 2014J the NIR light curve can yield an inde¬ 
pendent check on the nickel mass and a direct comparison to 
the direct measurements of the 7 —rays from the nickel decay 
(Churazov et al. 2014; Diehl et al. 2014). A calibration of a 
fair set of (preferably bolometric) peak luminosities and the 
derived nickel masses with our parameter t 2 should lead to 
the corresponding relation. 

We attempted to improve the IR first maximum for dis¬ 
tance measurements and looked for a possible correlation of 
Ml with the phase of the second maximum. There is a slight 
improvement in the J band and none in Y and H. The re¬ 
sulting scatter after correction in J band is 0.16 mag. This 
is comparable to the observed value from previous studies 
(Mandel et al. 2009; Kattner et al. 2012; Weyant et al. 2014). 
It remains to be seen whether larger samples will provide a 
better handle in the future. 

Since we specihcally studied the second maximum in the 
NIR light curves, we excluded SNe la, which do not display 
this feature. These are objects similar to either SN 1991bg, 
SN 2000cx or SN 2002cx and are typically faint and peculiar 
SNe (Filippenko et al. 1992; Leibundgut et al. 1993; Li et al. 
2001, 2003; Foley et al. 2013). They appear to also display 
a fainter hrst maximum in the NIR (Krisciunas et al. 2009). 
It would be interesting to check whether the uniform decline 
rate observed in the near-IR also applies to SNe la which do 
not show the second maximum. This could be used to check 
whether these explosions share some physical characteristics 
with the ones discussed here. There are only very few SNe la 
of this type and we do not find a conclusive answer. Increas¬ 
ing the number of SNe la with this information is important 
to assess the physical differences among the different SN la 
groups. The recent CFAIR2 catalogue (Friedman et al. 2014) 
will provide some of these data. 

The NIR light curves display a decline rate after the 
second maximum, which is significantly larger than the op¬ 
tical light curves at the same phase (e.g. Leibundgut 2000). 
At very late phases (^ 200 days) the near-IR light curves 
become nearly fiat as observed for SN 2001el, while there is 
no observable change in the optical decline rates (Stritzinger 
& Sollerman 2007). It would be interesting to observe the 
change in decline rate between 100 and 300 days. Pre¬ 
sumably, the internal structure of the explosions sets the 
transition towards the positron decays as the dominant en¬ 
ergy source, when the ejecta have thinned enough that the 
7 —rays escape freely. This phase might be correlated with 
other physical parameters, like the nickel and ejecta mass, 
determined through early light curves. 

A possible extension of this photometric study with 
detailed spectroscopic observations and theoretical spectral 
synthesis calculations might be worthwhile to check on the 
emergence of the various emission lines, trace the exact tran¬ 
sition to the flatter IR light curves and determine whether 
it indicates any differences in the structure of the SNe, e.g. 
transition to positron channel. 

Finally, Euclid will discover many SNe at near-NIR 
wavelengths out to cosmologically interesting redshifts (e.g. 
Astier et al. 2014). With the small scatter of the peak lumi¬ 


nosity, these observations will provide distances with largely 
reduced uncertainties due to reddening. Our study confirms 
the promise the NIR observations of SNe la offer. 
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Table Al. Y filter apparent Magnitude and phase of the two maxima and the minimum for the complete sample. The phases are given 
in days relative to the B maximum. . The full table is available online 


SN 

tl 

mi to 

mo 

t2 

m2 


(days) 

(days) 


(days) 


2004eo 

2004ey 

-4.48 ± 0.82 

15.90 ±0.03 


31.77 ± 0.20 

15.51 ± 0.01 

2004gs 


8.10 ±0.76 

17.69 ±0.02 

22.35 ±0.07 

17.16 ±0.03 

2004gu 




37.04 ±0.13 

17.67 ± 0.03 

2005A 


13.70 ±1.61 

17.69 ±0.41 

27.39 ±0.07 

16.00 ±0.11 


Table A2. J filter apparent magnitude and phase of the two maxima and the minimum for the complete sample. The phases are given 
in days relative to the B maximum. The full table is available online 


SN t\ mi to ^0 ^2 

(days) (days) (days) 


m2 


1980N 

1981B 

1986G 

1998bu 

1999ac 


21.86 ±2.70 
32.19 ±0.10 
16.40 ±1.40 
29.84 ±0.40 

27.00 ±2.00 15.61 ±0.50 


Table A3. H filter apparent magnitudes and phases of the two maxima and the minimum for the complete sample. The phases are 
given in days relative to the B maximum. The full table is available online 


SN ti mi to mo t2 

(days) (days) (days) 


m2 


1981B 

1986G 

1998bu 

1999ee 

2000E 


31.98 ±0.50 
12.70 ± 1.30 
29.18 ± 1.65 
30.50 ± 0.70 
25.10 ± 0.70 


10.01 ± 0.10 

15.11 ± 0.03 
13.85 ±0.10 


Table A4. Late time decline rate and luminosity at ±55days in Y JH. The full table is available online 


SN 

slope 

(mag/day) 

err 

M|55 

err 

Y filter 

2005M 

0.054 

0.002 

-17.73 

0.08 

2005el 

0.056 

0.002 

-16.72 

0.11 

2005na 

0.052 

0.003 

-17.57 

0.13 

2006D 

0.050 

0.001 

-16.62 

0.17 

2006X 

0.054 

0.003 

-16.67 

0.10 


Table A5. t 2 and epoch at which SN enters Lira law. The full table is available online 


SN 

^2 

(days) 

err 

(days) 

err 

Y filter 
2004gs 

22.34 

0.07 

24.17 

3.21 

2005am 

21.73 

0.12 

22.69 

1.26 

2005el 

24.96 

0.11 

25.78 

1.00 

2005na 

31.77 

3.46 

32.69 

1.26 

2006D 

24.88 

0.02 

25.00 

0.26 
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